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Abstract 

Using a linear mass spectrum of an SU(4) meson hexadecuplet, we derive 
the Gell-Mann-Okubo mass formula for the charmed mesons, in good agree- 
ment with experiment. Possible generalization of this method to a higher sym- 
metry group is briefly discussed. 
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The hadronic mass spectrum is an essential ingredient in theoretical investigations of 
the physics of strong interactions. It is well known that the correct thermodynamic de- 
scription of hot hadronic matter requires consideration of higher mass excited states, 
the resonances, whose contribution becomes essential at temperatures ~ 0(100 MeV) 
H], . The method for taking into account these resonances was suggested by Belenky 
and Landau || as considering unstable particles on an equal footing with the stable 
ones in the thermodynamic quantities; e.g., the formulas for the pressure and energy 
density in a resonance gas readQ 



_ _ mjT 2 ^ frnA 



pi, Pi = T—- Pil (2) 

where are the corresponding degeneracies (J and / are spin and isospin, respec- 
tively), 



7T 4 



Qi = — x 
y 90 



(2Jj + l)(2Jj + 1) for non — strange mesons 

4(2Jj + 1) for strange (K) mesons 

2(2 Ji + 1)(2I« + 1) x 7/8 for baryons 



These expressions may be rewritten with the help of a resonance spectrum, 

[■m.2 m 2 T 2 fm\ 

V = J m dm T(m)p(m), p(m) = ~7^r K 2 yfj > ( 3 ) 

ma . . ^dpim) . . . . 

dmr(m)p{m), p(m) =T P\ m )-> [V 

mi dl 



normalized as 



J mi 



dm r(m) =J2du (5) 



where m\ and mi are the masses of the lightest and heaviest species, respectively, 
entering the formulas (1),(2). 



1 For simplicity, we neglect the chemical potential and approximate the particle statistics by the 
Maxwell-Boltzmann one. 
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In both the statistical bootstrap model [||, Q and the dual resonance model |J , a 
resonance spectrum takes on the form 



r(m) ~ m a e m/T °, (6) 

where a and To are constants. The treatment of a hadronic resonance gas by means 
of the spectrum (6) leads to a singularity in the thermodynamic functions at T = T 
||, [5| and, in particular, to an infinite number of the effective degrees of freedom in 
the hadron phase, thus hindering a transition to the quark-gluon phase. Moreover, 
as shown by Fowler and Weiner J7[, an exponential mass spectrum of the form (6) 
is incompatible with the existence of the quark-gluon phase: in order that a phase 
transition from the hadron phase to the quark-gluon phase be possible, the hadronic 
spectrum cannot grow with m faster than a power. 

In our previous work || we considered a model for a transition from a phase 
of strongly interacting hadron constituents, described by a manifestly covariant rel- 
ativistic statistical mechanics which turned out to be a reliable framework in the 
description of realistic physical systems 0, to the hadron phase described by a res- 
onance spectrum, Eqs. (3), (4). An example of such a transition is what might be 
considered to be a relativistic high temperature Bose-Einstein condensation studied 
by the authors in ref. ||10||, which corresponds, in the way suggested by Haber and 



Weldon [II], to spontaneous flavor symmetry breakdown, SU(3)f — > SU(2)j x [7(l)y, 
upon which hadronic multiplets are formed, with the masses obeying the Gell-Mann- 
Okubo formulas |L2] 

- Y 2 



m = a + bY + c 



(7) 



here / and Y are the isospin and hypercharge, respectively, i is 2 for mesons and 1 
for baryons, and a, b, c are independent of I and Y but, in general, depend on (p, q), 
where (p, q) is any irreducible representation of SU(3). Then only the assumption on 
the overall degeneracy being conserved during the transition is required to lead to 
the unique form of a resonance spectrum in the hadron phase: 

r(m) = Cm, C = const. (8) 



Zhirov and Shuryak have found the same result on phenomenological grounds. 



As shown in ref. ||13| , the spectrum (8), used in the formulas (3), (4) (with the upper 
limit of integration infinity), leads to the equation of state p, p~T 6 ,p = p/5, called 
by Shuryak the "realistic" equation of state for hot hadronic matter 0, which has 
some experimental support. Zhirov and Shuryak [^] have calculated the velocity of 
sound, c 2 = dp/ dp = c 2 (T), withp and p defined in Eqs. (1),(2), and found that c 2 s {T) 
at first increases with T very quickly and then saturates at the value of c 2 s ~ 1/3 if 
only the pions are taken into account, and at c 2 s ~ 1/5 if resonances up to M ~ 1.7 
GeV are included. 
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We have checked the coincidence of the results given by the linear spectrum (8) 
with those obtained directly from Eq. (1) for the actual hadronic species with the cor- 
responding degeneracies, for all well-established multiplets, both mesonic and bary- 
onic, / 2 (1525), and found it excellent ||. Therefore, the theoretical conclusion that 
a linear spectrum is the actual spectrum in the description of individual hadronic 
multiplets finds its experimental confirmation as well. In our recent paper [14[] we 
applied a linear spectrum to the problem of establishing the correct qq assignment 
for the problematic meson nonets, like the scalar, axial-vector and tensor ones, and 
separating out non-gg mesons. 

The easiest way to see that a linear spectrum corresponds to the actual spectrum 
of a meson nonet is as follows^ Let us calculate the average mass squared for a spin-s 
nonet: 

, 2 \ _ YnOM 3m; + 4m; /2 + mjj, +mjj„ 
{mh = ^~- 9 ' (9) 

where mi, mi/2, mo, m ' are the masses of isovector, isospinor, and two isoscalar 
states, respectively, and the spin degeneracy, 2s + 1, cancels out. In general, the 
isoscalar states^ u > and u " are the octet u>s and singlet u mixed states because of 
SU(3) breaking, 

= u 8 cos 9m — sin 9m, 

uj " = u 8 sin 9m + w cos 9m, 

where 9m is a mixing angle. Assuming that the matrix element of the Hamiltonian 
between the states yields a mass squared, i.e., mi, = (ui '\H\u ') etc., one obtains 
from the above relations [|IJ|, 

ml =ml C os 2 9 M + ml S m 2 9 M -2ml 8 sm9 M cos9 M , (10) 

m 2 ,, = mj sin 2 9 M + m 2 cos 2 9 M + 2m 2 8 sin 9 M cos 9 M . (11) 
Since u ' and u » are orthogonal, one has further 

m o'o" = = ( m 8 — m o) sm ^M cos9m + mg 8 (cos 2 9m — sin 2 9 M )- (12) 

Eliminating mo and mos from (10)-(12) yields 

„ ml — m 2 , , 
tan 2 e M = ^ \. 13 

m Q " ~ m 8 

It also follows from (10), (11) that, independent of 9m, rn 2 , + m 2 „ = m\ + m\, and 
therefore, Eq. (9) may be rewritten as 

3m 2 + Am\ /2 + m\ + m 2 
(m )9 = '-g • (14) 

2 For a baryon multiplct, it is more difficult to show that the mass spectrum is linear, since the 
Gell-Mann-Okubo formulas are linear in mass for baryons. More detailed discussion is given in S. 



3 The u> ' is a mostly octet isoscalar. 
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For the octet, (3 mi, 4 mu2, 1 mg), the Gell-Mann-Okubo formula (as follows from 
(7)) is 

4mw 2 = 3jtij + mj. (15) 
Therefore, the average mass squared for the octet is 

3m 2 + 4m 2 /2 + m 2 , m 2 + m 2 
(™ )s = J} = 2 — ' ^ 

where Eq. (15) was used. In the exact SU(3) limit where the u, d and s quarks have 
equal masses, all the squared masses of the nonet states are equal as well. Since in 
this limit all the squared masses of the octet states are equal to the average mass 
squared of the octetQ, Eq. (16), the mass of the singlet should have the same value,Q 
i.e., 

Q m? + mo , . 

< = \ 8 - (17) 

With Eq. (15), it then follows from (17) that 

m + m l = 2m l/2i 

which reduces, through ??1q + m| = rriii + m^„ , to 

m 2 / + m^n = 2my 2 , (18) 
which is an extra Gell-Mann-Okubo mass relation for a nonet. We have checked this 



relation in a separate paper [IB], and found that with the experimentally available 
meson masses, the relative error in the values on the l.h.s. and r.h.s. of Eq. (18) does 
not exceed 3% for all well-established nonets (except for the pseudoscalar nonet for 

4 In a manifestly covariant theory, this holds since a total mass squared is rigorously conserved. 
In the standard framework, for pseudoscalar mesons, this is easily seen by using the lowest order 
relations [jl6|| mf = = 2mB, m/^, 2 = rn 2 K — (m + m s )B 1 where m — (m u + m<j)/2, and B is 
related to the quark condensate. Therefore, it follows from (15), (16) that m§ = 2/3 (2m s + m)B, 
(m 2 )g = 2/3 (m s + 2m)B = 2/3 (m u + rrid + m s )B. In the exact SU(3) limit, m u — nid = m s — fh, 
and hence m\ = m^ 2 = m| = (rn 2 )s = 2fhB. For higher mass mesons, since the states with equal 
isospin (and alternating parity) lie on linear Regge trajectories, one may expect the relations of 
the form (c = C/B) m\ — 2mB + C = (2m + c)B, m^ 2 = (m + m s )B + C = (m + m s + c)B, 

m\ = 2/3 (2m s +m)B + C = 2/3 (2m s +rn+3/2 c)B, consistent with the Gell-Mann-Okubo formula 
(15), leading to to 2 = m^ 2 — m| = (rn 2 )^ — 2mB + C in the SU(3) limit m u = md = m s = fh. 
For vector mesons, such a relation was obtained by Balazs in the flux-tube fragmentation approach 
to a low-mass hadronic spectrum ]T^ |, m 2 p — m 2 + l/2a , with a being a universal Regge slope, in 
good agreement with the experiment. 

5 It is also seen from the relations of a previous footnote: since the total mass squared of a nonet is 
proportional to the total mass of quarks the nonet members are made of, ^ i Qi m 2 = (12m+6m s )B+ 
9C, it follows from the above expressions for m 2 , m 2 ^ 2 and m| that m\ = 2/3 (2m + m s )B + C = 
(m 2 ) 8 = (m 2 ) 9 . 
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which Eq. (18) does not hold, perhaps because the r/o develops a large dynamical mass 
due to axial £7(1) symmetry breakdown before it mixes with the rjs to form the physical 
7] and r] states). For a singlet-octet mixing close to "ideal" one, tan 9 m ~ l/v2; it 
then follows from (13) that 

2mjy + itIqii ~ 3mg, 
which reduces, through (15), (18), to 

m » ~ mi. (19) 

Now it follows clearly that the ground states of all well-established nonets^ (except 
for the pseudoscalar one) are almost mass degenerate pairs, like (p, u;).[] In the close- 
to-ideal mixing case, Eq. (18) may be rewritten, with the help of (19), as 

ml + m\, ~ 2m\ /2 . (20) 

This relation for pseudoscalar and vector mesons with the ground states being the 
mass degenerate pairs (7r,^o) an d (p, u), respectively, was previously obtained by 
Balazs and Nicolescu using the dual-topological-unitarization approach to the con- 
finement region of hadronic physics (Eq. (21) of ref. [19|]). With (16) and (17), Eq. 
(10) finally reduces to 

(m 2 h = (21) 

which, of course, coincides with both, (m 2 )§ in (16) and m\ in (17), which is the 
property of the SU(3) limit (or the conservation of a total mass squared in a manifestly 
covariant theory). 

For the actual mass spectrum of the nonet, the average mass squared (9) may be 
represented in the form| 

. 9l f™ 8 dm rim) m 2 , s 

™ 9 = L 8 , ,\ , 22 

and one sees that the only choice for r(m) leading to the relation (21) is r(m) = Cm, 
C = const. Indeed, in this case 



f™* dm m 3 (ml — m\)/A m\ + m\ 
Q dm m = [mj - m\)/2 = 2 ' 



in agreement with (21). 

6 This is also true for qq assignment of the scalar meson nonet suggested by the authors in ref. 

s 

7 It follows from the relations of footnote 4 that, in the close-to- ideal mixing case, mjy ~ 2m s B+C 
and nitf, ~ 2mB + C = m\. 

8 Since m s > to, mi < mya < nig, as seen in the relations of footnote 4. Moreover, toi < too < 
TOg , and therefore, the range of integration in Eq. (22) is (mi, ms). 
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Evidently, one may choose an opposite way, viz., starting from a linear spectrum 
as the actual spectrum of a nonet, to derive the Gell-Mann-Okubo mass formula. To 
this end, one should first calculate the average mass squared, Eq. (21). Then one has 
to place 9 nonet states in the interval (mi, m 8 ) in a way that preserves the average 
mass squared. As we already know, the isoscalar singlet mass squared should coincide 
with the average mass squared; for the remaining 8 states one would have the relation 
(16) which would in turn reduce to the Gell-Mann-Okubo formula (15). One sees 
that the assumption of a linear mass spectrum turns out to be a good alternative to 
group theory assumptions on the form of the mass splitting, for the derivation of the 
Gell-Mann-Okubo type relations, which may be rather difficult technical task for a 
higher symmetry group. 

We now generalize the above derivation of the Gell-Mann-Okubo mass formula to 
the case of 4 flavors, by simply adding one more quark (c-quark). Then, in addition 
to the 9 states of a nonet, we have to add 7 more states, to finally form an SU(4) 
hexadecuplet: 4 cu,uc,cd,dc, 2 cs,sc, and 1 cc. Suppose that SU(3) symmetry is 
exact but the underlying SU(4:) one is broken by the c-quark mass. Then the 9 nonet 
states have equal mass squared which coincide with Wq, and of the remaining 7 states, 
6 (which are the combinations of c-quark with one of u, d, s) have equal masses as 
well. Now we have to place 7 states in the interval (mo, m c5 ) in a way that preserves 
the average mass squared (q stands for one of u, d, s) : 

6m^ + m^ _ m 2 + m 2 c5 _ 
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therefore 



(23) 



I2m% = 5m 2 c - c + 7ml (24) 

which is the generalization of the Gell-Mann-Okubo mass formula to the case of 
broken S77(4) but exact SU{3) symmetry. 

In a real world, both, S77(4) and SU(3), are broken, so that we have to distinguish 
between the masses of cu, cd, and cs states. Introducing the standard notations, 
D(cu,cd), D s (cs), \I/(cc), we have to modify Eq. (23), as follows: 

4m^ + 2m 2 Da + ml _ mg + m| 

7 _ 2 ; [Zb) 

thus 

8m 2 D + Am 2 Da = 5m| + 7m 2 0} (26) 

which is the Gell-Mann-Okubo mass formula for an SU(4) hexadecuplet (which has 
to be accompanied by those for an SU(3) nonet, Eqs. (15), (16)). Let us check this 
formula for the experimentally established masses of the charmed mesons, for the 
following multiplets (which are the only multiplets for which these mesons have been 
discovered experiment ally) []: 

9 The values of m\ given below are calculated from Eq. (17), which reduces, through (15), to 
ml = (mf + 2m 2 1/2 )/3. 
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1) 1 J PC = 0"+ pseudoscalar mesons, m(D) = 1.87 GeV, m{D s ) = 1.97 GeV, 
m(j] c ) = 2.98 GeV, and m\ = 0.17 GeV 2 . Therefore, one obtains 43.5 GeV 2 on the 
l.h.s. of Eq. (26), vs. 45.6 GeV 2 on the r.h.s. 

2) 1 3 Si J PC = 1— vector mesons, m(D*) = 2.01 GeV, m(D*) = 2.11 GeV, 
m(J/-9) = 3.09 GeV, and m 2 = 0.73 GeV 2 . Therefore, one has 50.1 GeV 2 on the 
l.h.s. of Eq. (26) vs. 52.8 GeV 2 on the r.h.s. 

3) 1 1 P 1 J PC = l+~ pseudovector mesons, m{D x ) = 2.42 GeV, m{D sl ) = 2.54 
GeV, m(h c (lP)) = 3.53 GeV, and m\ = 1.58 GeV 2 . In this case one has 72.7 GeV 2 
on the l.h.s of (26) vs. 73.3 GeV 2 on the r.h.s. 

4) 1 3 P 2 J PC = 2++ tensor mesons^, m(D* 2 ) = 2.46 GeV, m(D* 2 ) = 2.57 GeV, 
m( X c2(lP)) = 3.55 GeV, and m 2 = 1.94 GeV 2 . Now one has 74.8 GeV 2 on the l.h.s 
of (26) vs. 76.5 GeV 2 on the r.h.s. 

Thus, in either case, the relative error does not exceed 5%, in the third case the 
results almost coincide. 

One may try to generalize the presented method of the derivation of the Gell- 
Mann-Okubo formula from a linear mass spectrum to a higher symmetry group. Let 
us briefly discuss this point. Suppose, as previously, that SU(N) flavor symmetry is 
exact but the underlying SU(N + 1) one is broken by the mass of the (N + l)-quark. 
Then we have a mass degenerate SU(N) multiplet (N 2 states), and 2N + 1 more 
states, of which 2N are mass degenerate as well. Now we have to place these 2N + 1 
states in the mass interval (m(jv,jv), m (N+i,N+i)) in a way that preserves the average 
mass squared: 

2Nm 2 N ^ N+1 - ) + m 2 JV+ljJV+1 ) ^ rn 2 {N,N) + m (jV+i,iv+i) (C)7 , 
2N+ 1 ~ 2 ' 

This is the analog of the Gell-Mann-Okubo mass formula for an SU(N + 1) multiplet. 
It one writes down roughly (m^v and irtN+i are the masses of N- and (iV + l)-quarks, 
respectively) 



(n,n) - 2m N , m(Ar i7V+ i) ~ m N + m N+1 , m( N+ljN+1) ~ 2m N+1 (28) 



and neglects m% in comparison with m% +1 , one obtains 

2N 

m N+ x ~ — _ ^ m N . (29) 

Now it becomes clear why the quark masses are exactly what they are. If one goes 
from SU(3) to SU(4), one obtains from (29) ~ 3m3, i.e., m c ~ 3m s , where m s is 
the constituent s-quark mass. With m s ~ 0.5 — 0.55 GeV, one has m c ~ 1.5 — 1.6 



10 Although the Z?* 2 (2573) meson was omitted from the recent Meson Summary Table as "needs 
confirmation" , it is the best candidate on the 1 Z P^ 2 ++ cs-state since its width and decay modes 
are consistent with J p = 2 + |20|. 
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GeV, exactly the value which may be expected from the naive quark model and which 
is indicated by the recent Particle Data Group f2l| . If one further goes from SU(4) to 
SU(5), one obtains m 5 ~ 8/3 m 4 , i.e., m b ~ 8/3 m c . With m c ~ 1.5 — 1.6 GeV, this 
implies m b ~ 4.0 — 4.3 GeV, again in agreement with the value 4.1 — 4.5 GeV given by 
the recent Particle Data Group ^TJ and the expectations from the naive quark model. 



We therefore expect the corresponding Gell-Mann-Okubo formula for an SU (5) 25- 
plet to be in a fair agreement with the experiment. Unfortunately, if one goes from 
SU(5) to SU(6), one gets the value of the top-quark mass m t ~ 5/2 m& ~ 10 — 11 
GeV, which is one order of magnitude less than that expected from the Standard 
Model[^] ||21|| . The explanation of this fact may lie in the mass spectrum deviating 
from a linear form for the masses of the order of ~ 0{m b ). (E.g., the mass spectrum 
may be of the form r(m) ~ m ■ f(m), where f(m) is close to unity for the masses up 
to ~ 0(mb) and begins to grow rapidly for higher masses.) In fact, in a manifestly 
covariant theory, the mass spectrum is r(m) ~ m ■ exp(am 2 ), but normally one 
expects a to be very small. We also note that the actual mass spectrum cannot 
be of the Hagedorn form (6) because the exponent a in Eq. (6) is always negative 
(and is related to the number of transverse dimensions of a string theory ||23||), and 



therefore cannot be equal to unity. More detailed discussion on this point, as well 
as the Gell-Mann-Okubo formulas for SU(4) baryon multiplets, will be given in a 
separate publication. 
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